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SUMMARY 

Numerical simulations of the 3D MHD-equations that describe rotating magnetocon- 
vection in a Cartesian box have been performed using the code NIRVANA. The char- 
acteristics of averaged quantities like the turbulence intensity and the turbulent heat 
flux that are caused by the combined action of the small-scale fluctuations are com- 
puted. The correlation length of the turbulence significantly depends on the strength 
and orientation of the magnetic field and the anisotropic behavior of the turbulence 
intensity induced by Coriolis and Lorentz force is considerably more pronounced for 
faster rotation. The development of isotropic behavior on the small scales - as it is 
observed in pure rotating convection - vanishes even for a weak magnetic field which 
results in a turbulent flow that is dominated by the vertical component. In the presence 
of a horizontal magnetic field the vertical turbulent heat fiux slightly increases with 
increasing field strength, so that cooling of the rotating system is facilitated. Horizontal 
transport of heat is always directed westwards and towards the poles. The latter might 
be a source of a large-scale meridional fiow whereas the first would be important in 
global simulations in case of non-axisymmetric boundary conditions for the heat flux. 

Key words: magnetoconvection - anisotropic turbulence - geodynamo - turbulent 
heat fiux - MHD-simulations 



1 INTRODUCTION 

Anisotropic turbulence is a fundamental feature in many 
astro- or geophysical systems. Well known realizations are 
the turbulent motions in the solar convection zone or the 
flow of liquid iron in the outer core of the Earth. Both kinds 
of convectively driven flows are assumed to be responsible 
for dynamo action. Flow and field are attuned in a complex 
non-linear system where the flow produces the magnetic fleld 
which in turn backreacts on the field producing flow. Di- 
rect observations of the sun ofi'er the possibility to examine 
the details of the complicated interactions between convec- 
tion and magnetic flelds. On the top of the solar convection 
zone, in the photosphere, where rotational effects are rather 
unimportant, occasionally strong localized and radial ori- 
ented magnetic fields inhibit the convective transport of en- 
ergy which results in a cooler (and therefore darker) area, a 
sunspot (see e.g. Weiss, 1990). Numerical simulations of non- 
linear magnetoconvection confirm that a vertical magnetic 
fleld suppresses turbulent motions as well as it considerably 
reduces the size of the convection cells, a result that has been 
known from the linear stability analysis of Chandrasekhar 
(1961). The smaller cells come along with reduced varia- 
tions of the temperature so that the correlations of velocity 
and temperature fluctuations, the turbulent heat flux, de- 



crease with increasing field (Cattaneo et al., 2003). Towards 
the edge of the sunspot the field changes its orientation to 
a more horizontal direction. As a result the flow pattern 
changes to small brighter and darker elongated filaments 
that surround the sunspot (Weiss et al., 2004). Thermody- 
namic properties like the heat transport and the tempera- 
ture distribution are undoubtedly linked to the anisotropy 
induced by the dominant magnetic fleld component. Closely 
related to the interaction of convection, magnetic fleld and 
turbulent heat transport are questions about the existence 
of warmer (and therefore brighter) rings around sunspots 
(Eschrich & Krause, 1977; Riidiger & Kitchatinov, 2000) or 
the behavior of (star-) spots in fast rotating stars. In con- 
trast to the sun, observations of fast rotating stars using 
Doppler-Imaging techniques reveal an accumulation of star 
spots close to the poles (Strassmeier, 2002, 2006). This phe- 
nomenon clearly indicates the diflerent behavior of convec- 
tively driven turbulence in fast rotating objects and points 
out that a different dynamo mechanism is running where 
differential rotation is of minor importance (Bushby, 2003). 
Comparable considerations might also be of interest for the 
Earth because observations of the radial fleld component at 
the core mantle boundary as well as highly resolved simula- 
tions of the geodynamo show strong localized flux patches 
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whose pattern and pair-like occurrence resemble the behav- 
ior of sunspots (Roberts & Glatzmaier, 2000; Jackson, 2003) . 

However, the particular properties of a fast rotating 
planetary body like the Earth lead to physical conditions 
that differ significantly from the sun or stellar objects. A 
small density contrast between the top and the bottom of 
the fluid outer core and a nearly incompressible fluid result 
in a small Mach number flow which is essentially influenced 
by the non-linear back-reaction of a strong magnetic field. 
Convection in the fluid outer core is mainly affected by two 
dominant forces, the Lorentz force due to the strong mag- 
netic field and the Coriolis force due to the fast rotation of 
the Earth. In such a rapidly rotating convection driven dy- 
namo the turbulence is subject to three preferred directions 
defined by magnetic field, rotation axis and gravity. The 
arising anisotropics result in a plate like form of convection 
cells which are aligned along the rotation axis and elongated 
in the direction of the dominant magnetic field component 
(Braginsky & Meytlis, 1990; St. Pierre, 1996; Matsushima 
et al., 1999). 

Despite the progress in computational and numerical 
techniques that have been made during the recent years, it 
is still impossible to simulate turbulent convection in rotat- 
ing spherical shells with a sufficient resolution so that the 
resolved scale range of the turbulence is rather restricted 
(HoUerbach, 2003). Furthermore, anisotropy usually is ne- 
glected by an implicit assumption of a large eddy simulation 
where scalar parameters resemble the turbulent values of the 
diffusivities. Beside the fact that mostly the diffusivities - for 
the reason of numerical stability - have to be chosen much 
larger than even the (poorly known) "real" turbulent val- 
ues, this oversimplifying assumption describes an isotropic 
transport of flux and an isotropic dissipation of energy. This 
is only justified if either the resolution is high enough so 
that all dynamical important modes are resolved or if no 
preferred direction originated by external forces, stratifica- 
tion or boundary conditions exists. Both, in general, are not 
the case. In order to include anisotropic effects in a more 
sophisticated way, turbulence models have to be considered 
where tensorial expressions for the diffusivities take care of 
the directional dependence of the turbulence and the influ- 
ence of the unresolved scales is parameterized in terms of 
the resolved large-scale fields. 

Tensor coefficients for the viscous and thermal diffusiv- 
ity of a fast rotating and a strong field model of the Earth's 
core have been derived by Phillips & Ivers (2001, 2003). 
They present expressions that describe enhanced diffusion 
along a dominant azimuthal magnetic field. Non-diagonal 
elements are neglected so that diffusion in the horizontal di- 
rections due to a coupling between vertical and horizontal 
components of the turbulent flow is not included. The results 
are intended for the use in pseudo-spectral codes but at this 
stage no applications in numerical investigations are avail- 
able so that the consequences for global simulations remain 
unknown. 

A different approach is examined by Buffett (2003) 
and Matsushima (2004, 2005). They present a subgridscale 
(SGS) model which is essentially based on the self- similarity 
of the turbulence. By comparison of a highly resolved direct 
numerical simulation (DNS) with a large eddy simulation 
based on the SGS model, they show that the SGS model is 
able to predict appropriate anisotropic heat- and momentum 



fluxes. An advanced version of the SGS model that includes 
the Lorentz- and induction terms is introduced by Matsui & 
Buffett (2005) . They apply the so called non-linear gradient 
model, an adaption of the similarity model that addition- 
ally is based on the local character of the turbulence. In 
general, the SGS models coincide rather well with the di- 
rect numerical simulations, however, both models are not 
able to reproduce all details of the spatiotemporal behavior. 
Furthermore, systematic deviations are obtained in mean 
quantities, like kinetic and magnetic energy, which are sit- 
uated between the values from the resolved DNS and the 
values of a truncated, unresolved DNS. 

In order to evaluate the functionality of SGS models it 
is useful to further investigate the development of a convec- 
tively driven turbulence in view of the dependence on field 
strength and latitude. Subject matter of the present paper is 
the connection between turbulence characteristics, that are 
caused by the average action of the small-scale fluctuations 
and the turbulent transport of heat which may essentially 
contribute to the conditions that determine the temperature 
distribution within the fluid outer core. A simplified Carte- 
sian system is examined where a convectively driven turbu- 
lence in a conducting fluid is subject to fast rotation and an 
externally imposed magnetic field. As we do not intend to 
reproduce geophysical realities in all details, the presented 
investigation is based on a highly idealized system where the 
emphasis lies on the anisotropy inducing effects of rotation 
and magnetic field. Properties like the exact behavior of the 
temperature gradient, the equation of state or the influence 
of curvature are assumed to be of secondary importance. 



2 THE MODEL 

2.1 General properties 

A detailed description of the applied local model can be 
found in Giesecke et al. (2005). Figure [1] shows a sketch of 
the computational domain, a Cartesian box placed some- 
where on a spherical shell at a co-latitude 6. The unit vec- 
tors x,y,z form a right-handed co-rotating Cartesian co- 
ordinate system with x pointing towards the equator, y 
pointing in the toroidal direction (from west to east) and 
z pointing from the bottom to the top of the box. In global 
spherical coordinates, z represents the radial direction r 
(oriented from inside to outside), y the azimuthal direc- 
tion (fi (oriented eastwards) and x the meridional direction 
6 (oriented towards the equator). The angular velocity 
in the co-rotating local box coordinate system is given by 
n = —SlsmOx + QcosOz where Q is the angular veloc- 
ity of the rotating spherical shell. The box with an aspect 
ratio 8:8:1 is placed at different latitudes on the northern 
hemisphere of the rotating spherical shell and a standard 
resolution of 100 x 100 x 80 grid points is used in all calcu- 
lations. The co-latitude angle 6 is varied from 6 = Q° (north 
pole) to e = 75°. 



2.2 Equations 

The MHD-equations for a rotating fluid, including the effects 
of thermal conduction, compressibility, viscous friction and 
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Figure 1. Model box as part of a rotating spherical shell at 
co-latitude 9. 



losses due to magnetic diffusivity, are 

dtp = -V-(pn), (1) 
dt{pu) = -V ■ {puu) -Vp + V ■ a + pg (2) 

+ — (V X B) X B- 2pS7 X u, 
Mo 

dte = -V-{eu)-pV-u+a-Vu + — \V X Bf (3) 

po 

+v-(Kvr), 

dtB = V X (u X B - r;V X B). (4) 

Here, p denotes the density, u the velocity, p the pressure, B 
the magnetic flux density, T the temperature and e the ther- 
mal energy density. A constant gravitational field g — ~gz 
is assumed within the domain. The viscous stress tensor g 
is given by aij = vp{djUi -\- diUj — 2/3V • u5ij). u denotes 
the kinematic viscosity and k, the thermal conductivity co- 
efficient. The values of k, the dynamic viscosity Vdyu = fp 
and the magnetic diffusivity 77 are constant within the box 
volume. An ideal gas equation of state is assumed: 



(7-l)e 



Jl 

mp' 



PT 



(5) 



where k is the Boltzmann constant, m the atomic mass unit, 
p, the mean molecular weight {p, = 1 for all runs) and 7 = 
Cp/cv = 5/3 is the ratio of Cp, the specific heat at constant 
pressure, and Cv, the specific heat at constant volume. The 
permeability po is given by the vacuum value po = 4:Tt x 
10-^VsA-im~\ 

The equations (1) - (4) are solved applying the code 
NIRVANA (Ziegler 1998, 1999). NIRVANA makes use of 
a dimensional- and operator splitting approach. A second- 
order accurate finite-volume scheme with a piecewise linear 
reconstruction and monotonized slope limiter (van Leer) is 
used for the hydrodynamic advection part of the equations. 
A constraint transport solver utilizing the method of charac- 
teristics is employed for the numerical solution of the induc- 
tion equation. The time integration of the source terms is 



performed by an explicit Euler scheme (except the Coriolis 
force which is treated analytically). 



2.3 Initial state and input parameters 

The simulations are started with an initial state that is de- 
termined by a hydrostatic equilibrium dzP = —pg and a 
polytropic temperature distribution T — Tq (p/po)^ in the 
absence of motions (the subscript refers to values taken at 
the top boundary of the box). With the equation of state 
the initial density distribution is given by 



Po 



i/r 



(6) 



where d stands for the vertical box extension and the poly- 
tropic index F is given by T = In (l -I- /To) /In^. The 
stratification index ^ — pbot/po, the temperature To and 
the global temperature gradient dT/dz are prescribed input 
parameters. 

The parameters v, k. and 77 are calculated from the 
Rayleigh number Ra, 



Ra = 



Cpd* pg f dT 



9 



T V dz 



(7) 



with Cp = k{mp)~^'y{'y — 1)~^ the specific heat at constant 
pressure, the Prandtl number Pr = vpCp/n and the magnetic 
Prandtl number Pm = v j-r). The basic parameter set used 
for all simulations that are presented in this paper is given 
by Ra = 10^, Pr = 0.5 and Pm = 0.5. The diffusivity param- 
eters are scalar quantities assuming that in the local model 
the resolved scales contain most of the energy, so that the 
anisotropic transport is described with sufficient accuracy 
and the remaining non-resolved modes are less important. 
The rotation rate ft is parameterized by the Taylor number 
given by 



Ta : 



(8) 



which is related to the Ekman number by Ta — Ek^ . Es- 
sentially, a rotating system is examined where Ta is set to 
a value of 10^. A few simulations with a Taylor number 
Ta = 10^ have been performed which demonstrate the sig- 
nificant changes in the behavior of the turbulence as the 
rotation rate increases. 

The magnetic field is expressed by the Elsasser number: 



A = 



2Qpoprj 



(9) 



A represents the relation of the Lorentz force to the Cori- 
olis force, which are assumed to be of the same order of 
magnitude within the fluid outer core. The original picture 
favors A ~ 0(1) because at this value the critical Rayleigh 
number -Rflcrit at which the onset of thermal convection oc- 
curs, becomes minimal. In that case the Coriolis force can 
be balanced by the Lorentz force. Then, there remains no 
need for balancing these terms with the (extremely small) 
viscous terms which would force the convection to occur on 
very short length scales. Such a flow would be much more 
difficult to maintain which is evident on the basis of the 
suppression of convective motions in simple rotating con- 
vection or non-rotating magnetoconvection. For a detailed 
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description of tlie underlying mechanism see e.g. Riidiger & 
Hollerbach (2004). However, Zhang & Jones (1994) did not 
find such a minimum for convection in a rotating spherical 
shell, but they were able to show that the existence of a sta- 
ble magnetic field requires an Elsasser number in the range 
l^A^ 10. 

For the comparison of the magnetic quenching character 
in simulations with different rotation rates, A is not a suit- 
able parameterization of the field strength. For that case, 
the dependence on the imposed field strength is expressed 
in units of the equipartition field strength B^q defined by 

-Bcq = y/JIopUruis, (10) 

where Urms denotes the root-mean-square velocity (see def- 
inition in the following section) for non-rotating and non- 
magnetic convection (in code units: Urms ~ 19.5). 

2.4 Averaging procedure 

Due to the small stratification all considered quantities do 
only weakly depend on the vertical coordinate z (except 
close to the upper and the lower boundaries). Therefore it is 
justified to characterize the turbulence properties of a fluc- 
tuating quantity by a volume average given by 

(^') = ^EC^>'^'^-W-7.w). (11) 

On the left hand side, /' represents a fluctuating quantity 
and the brackets {•) denote the averaging procedure. On 
the right hand side, /i,j,fe(t) represents the numerically com- 
puted value of the quantity / at a certain grid cell labeled 
by i,j,k at a certain time t. f^{t) represents the horizon- 
tal average of the considered quantity. A'' = nxnyUz denotes 
the number of grid cells used for averaging. Time averaging 
is done for periods with no significant changes in the sta- 
tistically steady convection state. Convergence of the time- 
averaged solutions has been checked by comparing results 
obtained from averages over different (increasing) periods. 
Reasonable results are obtained for averaging periods larger 
than at least 20 turnover times Tadv = d/urms, where the 
root mean square velocity Urms is defined as Urms = \/ {u''^). 

All quantities show fluctuations in time with a stan- 
dard deviation of the order of 10 percent of their mean 
value. These fluctuations become smaller for increasing field 
strength, but within the context of this paper they do not 
exhibit any essential features so that the standard deviations 
are omitted for the reason of clarity. 

2.5 Boundary Conditions 

All quantities are subject to periodic boundary conditions 
in the horizontal directions. At the top and at the bottom of 
the computational domain constant values for density and 
temperature are imposed. The vertical boundary condition 
for the magnetic field is a perfect conductor condition, and 
a stress-free boundary condition is adopted for the horizon- 
tal components of the velocity and Uy. Impermeable box 
walls at the top and the bottom lead to a vanishing at 
the vertical boundaries. Table [T] summarizes the boundary 



Table 1. Vertical boundary conditions 
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conditions and gives the initial values for density and tem- 
perature which describe the overall stratification and the 
global temperature gradient. For all simulations, tempera- 
ture and density at the top of the box are scaled to unity, as 
it is the case for the global temperature gradient dT/dz and 
the box height d. A stratification index of ^ = pbot/po = 1.1 
is used. 

The boundary conditions certainly influence the behav- 
ior of the flow within the computational domain. To reduce 
the undesired effects of the boundaries all volume averages 
are performed over the inner part of the computational do- 
main (between z — 0.2 to 0.8). Naturally, such a simplifying 
box model is not able to represent the large scale flows that 
may occur in a rotating sphere. Especially the existence of 
periodic boundary conditions in the horizontal directions fll- 
ters out the fraction of the solution with wavelengths larger 
than the horizontal extension of the box (of course, this be- 
havior is desired, because here only the small scale fluctua- 
tions are of interest). However, such filtering is not possible 
in the vertical direction due to the existence of density strat- 
ification and temperature gradient. 



3 RESULTS 

3.1 Convection pattern and anisotropy 

Within the applied parameter regime all simulations result 
in a pressure dominated (/? = 2fiop/B^ 3> 1) low Mach 
number fiow (Ma = Umis/cs 1, with the sound speed Cs = 
yJ^'Toj^—y Depending on rotation rate and field strength 
Ma is of the order of 10"^ (e.g. for A = 1 and Ta = 10^ the 
characteristic Mach number is given by Ma « 0.03). The 
corresponding Reynolds numbers Re = Urmsd/i^ are of the 
order of 100. The results that are presented in the following 
are obtained from simulations with an imposed horizontal 
field {By) which corresponds to a toroidal field in spherical 
coordinates. The typical structures of the convective mo- 
tions are visualized in Fig. (2] The upper row presents the 
vertical velocity pattern obtained from a simulation run with 
Ta = 10*^ (left: A « 1 respectively B = O.lBoq, right: A « 4, 
respectively B = 0.2i3cq). The lower row presents the case 
Ta = 10'^ (left: A « 1 respectively B = 0.2Boq, right: A « 4, 
respectively B = 0.4i3oq). AH plots show a time snapshot 
of the vertical velocity pattern from a simulation run at 
9 = 45°. 

The cell- like objects represented by the volume rendered 
iso-surfaces indicate coherent large-scale structures which 
span the full vertical extend of the domain. These convection 
cells are continously formed, rearranged and dissolved lead- 
ing to a quasi-stationary turbulent pattern which is more 
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Figure 2. z-componcnt of the velocity for Ta = 10'^ (upper row) and Ta = 10^ (lower row). Left side: A = 1, right side: A = 4. All 
plots present a snapshot at 6 = 45°. The solid iso-surfaces denote the vertical velocity at Uz = nrms, and the transparent iso-surfaces 
denote the pattern at = O.Snrms (upper left tirms ~ 11.5, upper right: u^ms ~ 10.8, lower left: nrms ~ 5.4, lower right: nrms ~ 4.7, 
all numbers in code units). A horizontal cut of the velocity field at a depth z = 0.5 is projected on a plane below the three dimensional 
box. Blue (dark) tones represent downwards oriented flow and yellow/red (light) colors represent upwards oriented flow. 



irregular in the interior of the domain. According to the 
Taylor-Proudman theorem the convection cells are aligned 
with the rotation axis. The average horizontal size of the 
cells strongly depends on the rotation rate as well as the 
magnetic field strength and direction. The cell size is re- 
duced for an increasing rotation rate and if the horizontal 
magnetic field exceeds a certain threshold the cells become 
elongated in direction of the dominant magnetic field com- 
ponent. With respect to Fig.[2]this is only noticeable in the 
faster rotating case whereas for Ta — 10^ the convection 



pattern remains nearly unaffected as the Elsasser number is 
increased from A ~ 1 to A ~ 4. 

Opposite to a configuration with a large density con- 
trast where the structure of the convection pattern consists 
of isolated, broad warm upfiows and narrow network-like, 
cold and strong downflows, within a weakly stratified layer 
the separation of up- and downflows is less pronounced. In- 
stead the upfiows form a kind of unconnected single plume- 
like structures enclosed by a loosely connected broad net- 
work of downflows. Up- and downflows are of approximately 



6 A. Giesecke 



equal amplitude and occupy roughly an equal area in the 
horizontal plane. 

To characterize the anisotropy of the turbulence and to 
quantify the visual impression obtained from Fig. [2] the av- 
erage extension of a convection cell is estimated in both hor- 
izontal directions. The two-point correlation function Q^z is 
defined for the direction perpendicular to the imposed mag- 
netic field by 



{u'^{x)u'^{x + 5x)) 

and for the direction parallel to the field: 

{u',{y)u',{y + 5y)) 



Qzz{5y) = 



{n'zivY) 



(12) 



(13) 



Qzz can be used as a convenient measure for the average 
horizontal extension of a convection cell. Motions within a 
cell are oriented in the same direction so that they are highly 
correlated whereas no correlations exist if the distance be- 
tween the considered two points (given by 5x, respectively 
Sy) exceeds the size of the cell. The mean cell size is inter- 
preted as the correlation length Acorr or Taylor microscale, 
which is the characteristic length-scale of the vorticity fila- 
ments observed in swirling fiows. Acorr is estimated from a 
function /, defined as 



f{Sx) = 1 - \J,,{5xf 



(14) 



which is adjusted to the decreasing part of the two-point 
correlation function Qzz- 

Figure O shows Qzz{5x) (upper panel) and Qzz{5y) 
(lower panel) for A = 0.1, 1, 10, 100 (solid curves) and the 
corresponding functions / (dashed curves). The shape of the 
curves of Qzz resembles the decrease of the correlation of the 
turbulent vertical velocity at two different coordinates with 
increasing distance between these two points. The behavior 
of Qzz{6x) which determines the correlation length perpen- 
dicular to By is independent of the imposed field strength 
(therefore the curves in the upper panel of Fig. [3] are not 
labeled by A). However, for large 5x the curves in the up- 
per panel exhibit an oscillating-like structure, which evolves 
in dependence of the field strength. With increasing field 
strength the amplitude increases and the length scale of 
these oscillations decreases. In case of weak fields the corre- 
lations outside of one cell are averaged out due to the dis- 
ordered occurrence of convection cells. For strong fields the 
sheetlike structures that can be observed on the lower right 
panel of Fig. 2 establish a more ordered structure which re- 
sembles in the oscillations of the correlation function in the 
upper panel of Fig. [S] 

The influence of the magnetic field on the cell struc- 
ture is considerably pronounced in the lower panel of Fig. [3] 
where Qzz is broadened with increasing fleld strength. The 
resulting correlation lengths A!^"" - estimated independently 
for both horizontal directions from Eq. p4p - are shown in 
Fig. |4l The black curve denotes the case Ta = 10^ where the 
solid (dotted) line represents Xi^'''' [\'^"). The gray curve 
shows the same quantities for Ta — 10^. As already indi- 
cated in the upper panel of Fig. [S] for both rotation rates, 
X^x^^ is nearly independent from the field strength so that 
the extension of the cell in x-direction is only determined 
by the rotation rate. According to the linear theory the pre- 
ferred length scale for the onset of rotating convection scales 




6y 

Figure 3. Time average of the two-point-correlation functions. 
The solid Unes represent Qzz{Sx), respectively Qzz{&y) and the 
dashed Unes represent a fit according to Eq. II14I I to the innermost 
part of Qzz- Ta = lO'^, = 0°, A « 0.1, 1, 10, 100. 
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Figure 4. Time and volume average of the correlation length 
A™J^'' in X- and j/-direction in dependence of the imposed magnetic 
field strength. The soUd (dotted) line denotes the extension in x- 
direction (j/-direction) for Ta = 10® (black) and for Ta = 10^ 
(gray). 6 = 0°. 
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as Ta"^'''' (for sufficient fast rotation). Scaling laws for rotat- 
ing finite amplitude convection have been examined by Stell- 
mach & Hansen (2004). They confirmed the above denoted 
scaling law for the size of a convection cell by evaluating the 
preferred wave number k = from the maximum value 

of the kinetic energy spectra. Within the rather restricted 
range of rotation rates our results reconfirm this scaling for 

Acorr 
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The correlation length parallel to the imposed field, 
Xy°'^'^ , resembles the increasing extension of the convection 
cells in direction of the imposed field. The field strength at 
which the transition to an anisotropic, elongated cell occurs 
only weakly depends on the rotation rate, however, the in- 
crease of A^""^"^ is slightly stronger for slower rotation. 

For increasing co-latitude 9 the convection cells become 
more and more tilted (according to the Taylor-Proudman 
theorem). For a higher co- latitude, A™'"^ increases due to 
this geometric effect, whereas A^""^"^ remains nearly unaf- 
fected (corresponding plots are omitted). 

The anisotropic character of the turbulence is also ap- 
parent in the behavior of the turbulence intensities {u'i^). 
For a quantitative examination the following functions are 
introduced: 
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(15) 



(see also Kapyla et al., 2004, and note the different defini- 
tion for Ay). The anisotropy of the turbulence intensities is 
described by Ah and Ay in the following way: 

{< dominated by turbulence perpendicular to 
= horizontal isotropy 
> dominated by turbulence parallel to By 

^ j < dominated by vertical turbulence 
^ I > dominated by horizontal turbulence 

As expected, in the presence of a magnetic field the mo- 
tions perpendicular By are more strongly suppressed than 
the motions parallel By. The resulting behavior of Ah at 
9 = 0° is shown in the upper panel of Fig. [5] which denotes 
the transition from the horizontal isotropic turbulence to 
an anisotropic state. A dependence on the rotation rate is 
evident. For Ta — 10® we approximately obtain horizontal 
isotropy up to field strengths corresponding to B ~ 0.4_Bcq 
(A « 10) whereas faster rotation results in An > except 
for very weak magnetic fields. 

The vertical anisotropy. Ay, is presented in the lower 
panel of Fig. \5\ For negligible magnetic fields an increase of 
the rotation velocity leads to an isotropic behavior on the 
small scales {Ay — > for _B — >■ 0), a property that has been 
found by Brummell et al. (1996). Already weak horizontal 
fields break this tendency and result in a turbulence which 
is dominated by the vertical component. Independent from 
the rotation rate the vertical anisotropy exhibits a maximum 
around By « O.SBoq- The dominance of the vertical compo- 
nent is characteristic for the behavior of the turbulence close 
to the pole and is not maintained for higher co-latitudes. The 
transition from vertical to horizontal dominated turbulence 
is determined by the orientation of the coherent cell struc- 
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Figure 5. Upper panel: Time and volume average of the horizon- 
tal anisotropy Ah for Ta = 10^ (solid line) and Ta = 10^ (dot- 
ted line). Lower Panel: Time and volume average of the vertical 
anisotropy Ay for Ta = 10® (solid line) and Ta = 10^ (dotted 
line). = 0°. 



tures, i.e. by the orientation along the rotation axis and by 
the elongation parallel to the dominant field component. 

3.2 Turbulent heat flux 

3.2.1 General description of the turbulent transport of 
heat 

For statistically steady convection the turbulent heat fiux 
Fr"^ is given by: 



5. 



pCp 



dT 
dxi 



9i 



(16) 



(see e.g. Riidiger, 1989) where mj denotes the tensorial heat 
conductivity which is related to the thermal diffusivity ten- 
sor Xij by Kij — pCpXij- In the presented model the large- 
scale temperature gradient dT/dxj and the gravity g are 
constant and oriented parallel to the z-axis so that we are 
only able to discuss the case j = z. Since, in addition, the 
density p remains approximately constant within the com- 
putational domain the simple relation Kiz oc {u'^T') holds to 
a good approximation and, in principal, immediately 
be determined from the numerically computed correlations 
{u[T'). It is convenient to define a normalized heat flux 



8 A. Giesecke 




Figure 6. Scatter plot of the normalized vertical heat flux in 
dependence of the turbulent up- and downflows at 2 = 0.5. Up- 
per left panel: convection, upper right panel: rotating convec- 
tion (Ta = 10^, 6 = 45°) , lower left panel: magnetoconvection 
{Bz ~ 0.2iJcq), lower right panel: rotating magnetoconvection, 
A = 1 (corresponds to By ^ 0.2Boq), Ta = lO'' and 6» = 45°. 



^ pconv / C^rji \ 

Fi = = {uiT = ^ (17) 

which is the quantity that is discussed with regard to the 
turbulent transport of heat in the remainder of this paper. 

The vertical heat flux is caused by turbulent upflows 
of warmer fluid and downflows of cooler fluid. Mainly the 
fastest up- and downflows contribute to the net flux which 
can be seen in Fig. [5] where time snapshots of the normal- 
ized flux u'^T' evaluated in a horizontal plane at z = 0.5 are 
presented in various scatter plots. Each dot denotes u'-^T' 
at a certain grid cell in dependence of the local vertical ve- 
locity u'.^. The sum of all values (divided by the number 
of data points) represents the horizontal average of {u'^T') 
at 2 = 0.5. The four panels show the results from simula- 
tions of convection (upper left panel), rotating convection 
(upper right panel), magnetoconvection (lower left panel) 
and rotating magnetoconvection (lower right panel). Qual- 
itatively the behavior is similar for all four rather different 
states which only distinguish in the amplitude of u'^T' and 
of the velocity u'^. The downflows (u^ < 0) and the upflows 
{u'z > 0) contribute with nearly the same amount to the net 
flux (although it seems that the downflows are slightly pre- 
vailing). This differs from the results for a strongly stratified 
layer presented by Ziegler (2002) who obtained a preponder- 
ance of the contributions by the downflows in non-rotating 
magnetoconvection which does not occur in sufficient fast 
rotating magnetoconvection. 

Typical volume averaged values of the vertical heat flux 
for moderate field strength and Ta = lO'^ are given by 
{u'^T') ~ 0.2 (in code units, see also section 3.2.2.). From 
Eq. (17) we therefore retrieve Xzz ~ 0-4, which is about 
four times larger than the value for the molecular diffusivity 
X = 0.1 as it results from the input parameters (all numbers 
in code units, see section 2.3 and 2.5). Although only moder- 
ate turbulence is examined, the obtained coefficient for the 
vertical transport of heat exceeds the molecular one, so that 



the transport of heat by turbulent motions is more effective 
than the molecular transport through heat conduction. 



3.2.2 Vertical heat flux and turbulence intensity 

From the quasi-linear theory a simple expression relating the 

thermal diffusivity tensor x s^nd Q^^ , the Fourier transform 

of the one-point correlation tensor Qij — {u'i(x,t)uj{x,t)) 

exists: 

[ /■ xfc^Q.,(fc,^) ^,^ 

For X ^ 0, the integrand can be replaced by a (5-function, 
so that: 

X^J =^ j Q^J{k,0)dk = i y Q^,{0,T)dr « ^TcorrQ^J (19) 

where the r- integral is approximated by Tcorr, the correla- 
tion time of the turbulence. A detailed explanation can be 
found in Riidiger (1989). This is a rather rough estimation 
which is based on the possibility to express the temperature 
fiuctuations T' through the fiuctuating velocity u' and is 
valid in the mixing-length approximation, where the turbu- 
lence spectra is dominated by one scale: ~ 5{k — t^I^)5{uj). 
Eq. (|19p might not be appropriate for a fast rotating system 
under the simultaneous influence of a magnetic field but 
it should deliver a general tendency for a relation between 
the turbulent velocity fluctuations and the (turbulent) ther- 
mal diffusivity tensor, respectively the normalized heat flux 

T-iconf 

In Fig. [7] the behavior of the vertical turbulence inten- 
sity {uf) (upper panel) is compared with the development 
of the normalized vertical heat flux, {u'^T') (lower panel) for 
Ta — 10^. Here, only the differences between the pole (solid 
line) and a co-latitude at 9 = 45° (dotted line) are pointed 
out. A more detailed analysis of the angular dependence is 
given in the next section (|3.2.4|l . As the most characteris- 
tic feature {u'^) increases with increasing field strength and 
exhibits a sharp drop for A ^ 2 at 61 = 0°. The maximum 
value for {u'^) is retrieved around A ~ 1...2. In this param- 
eter regime (where Lorentz force and Coriolis force are of 
the same order of magnitude) the turbulence is enhanced 
compared to rotating non-magnetic convection. This fea- 
ture does not appear at ^ = 45° or for a slower rotation 
rate (Fig. [8]). The enhanced turbulence around A ~ 1 has 
been predicted by Chandrasekhar (1961) who showed that 
the critical Rayleigh number Racrit is minimal if Lorentz and 
Coriolis forces are comparable. In case of a fixed (overcriti- 
cal) Rayleigh number a drop of Racrit leads to a more over- 
critical convection which results in a stronger driven flow and 
a more vigorous turbulence occurs. Although the linear sta- 
bility analysis of Chandrasekhar (1961) has been performed 
for a vertical oriented field (-Bz), the general trend might 
also be true for a horizontal field. From the linear analy- 
sis it is also known that this effect is more pronounced for 
faster rotating systems so that the maximum structure of 
(uf) will be more dominant for even faster rotation, a re- 
sult that already has been confirmed by Stellmach & Hansen 
(2004). However, the authors also obtained hints, that this 
effect - predicted from a linear stability analysis - vanishes 
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Figure 7. Time and volume average of the turbulence intensity 
and heat flux in dependence of the Elsasser number A. Upper 
panel: {uf), lower panel: {u'^T'). Ta = 10^. 



Figure 8. Time and volume average of the turbulence intensity 
{uf) and the turbulent heat flux (u'^T') for Ta = 10®. 



for stronger driven flows where the non-linearities dominate 
the final state. 

In accordance with the behavior of {u'^) the normal- 
ized vertical heat flux increases with the field strength be- 
low A ~ 1, but {u'^T') remains at a high level and only 
weakly decreases for strong fields. Although the fluctuating 
fluid motions are inhibited in the presence of a strong mag- 
netic field, this is obviously not the case for the temperature 
fluctuations, so that the correlation of these two quantities is 
not suppressed with increasing field strength. This feature is 
not observable at a higher co-latitude where {u'^T') exhibits 
a maximum around A « 0.5. 

A second distinctive feature is reflected in the crossover 
of the two curves for ^ = 0° and 6 = 45° in Fig. [T] For 
weak fields the turbulence intensity and the vertical heat 
flux a.t 6 = 45° are larger than at the pole whereas both 
quantities behave the other way round for a field strength 
above A « 0.1. ..0.5. Such a behavior does not occur for 
slower rotation (see Fig. [8]) or in case of a vertically oriented 
field (Fig. [9]). The vertical turbulence intensity and heat flux 
under influence of an external imposed Bz are illustrated 
in Fig. [9] (note, that in this case the magnetic boundary 
conditions have to be changed to = By = Q and dzB^ = 
0). Independent of the field strength, {u'^) and {u'^T') are 
slightly larger at 45° and both quantities show nearly the 
same dependence on the field strength. At the pole both 



quantities exhibit a weakly pronounced maximum around 
A ~ 1 and a very effective suppression occurs for A J> 2. 

3. 2. 3 Horizontal heat fluxes and Reynolds stresses 

In the following, again, a horizontal field is imposed and 
the horizontal transport of heat is discussed. The horizontal 
components of the normalized heat flux, {u'^T') and {u'yT') 
and the corresponding Reynolds stresses Qxz and Qyz van- 
ish at the poles. This expected behavior is nearly perfectly 
reflected in the numerical results and therefore the depen- 
dence on the field strength is only considered a.t 6 = 45°. 
Figure [To] shows Q^z and Qyz (upper panel) in comparison 
with (u'^T') and {u'yT') (lower panel). In all cases the consid- 
ered quantities are negative. Therefore, the latitudinal com- 
ponent, {u'^T'}, describes a transport of heat towards the 
(north-)pole. The absolute value increases with increasing 
field strength whereas the absolute value of Qxz decreases 
with increasing field strength. As a consequence the relation 
between {u'^T'} and Qxz cannot be described by a simple 
expression given by Eq. (fT^ . 

In a spherical geometry the y-component of the heat 
flux corresponds to an azimuthal transport of heat. {u'yT'), 
as it is also always negative, describes a westward directed 
transport of heat. The absolute value exhibits a distinct 
maximum value around A « 10 and for very strong fields 
{u'yT'} vanishes. The behavior of the azimuthal heat flux 
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roughly resembles the profile of Qyz, although there exists a 
small offset for the location of the minimum of Qyz towards 
weaker magnetic fields. 

3.2.4 Angular dependence 

In case of rotating turbulence without any other preferred 
direction the turbulent diffusivity tensor Xij is given by 

Xi] = XT<5ij + X| I ^^Jj^ + X^^P3 (20) 

(see e.g. Kitchatinov et al. 1994). x\\ describes enhanced 
thermal diffusion parallel to the rotation axis and the expres- 
sion x^ipj^p describes enhanced/reduced diffusion in the az- 
imuthal direction (which can be neglected in an axisymmet- 
ric configuration) . In the applied Cartesian system where we 
are restricted to a vertical large scale temperature gradient 
the following three elements of Xij can be computed: 

Xx^ = —XII sin 6* cos 6* (21) 
Xy^ = -X^? sin 61 (22) 
Xzz = XT+Xlicos^e* (23) 

These expressions neglect the influence of the magnetic field 
and might therefore be valid only for small values of A. 



The angular dependence of the Reynolds stresses and 
the normalized turbulent heat fiux are shown in Fig. [TT] (left 
hand side: Qiz; right hand side {u'iT')). The solid (dotted, 
dashed) lines denote the results for A = 0.1(1, 10) (corre- 
sponding to By = 0.07(0.2, 0.7)-Boq). In all cases the merid- 
ional component {u'^T') is negative. For A ~ 0.1 and A ~ 1 
the absolute value exhibits a maximum around 9 — 45° 
and in both cases {u'^T') can be described by a function 
oc sin 9 cos 9 as indicated by equation (|21|l . 

The azimuthal component {u'yT') is also always nega- 
tive and the absolute values increase towards the equator so 
that at least the tendency coincides with the behavior indi- 
cated by expression ((22}. \{u'yT')\ increases with increasing 
By which becomes more pronounced towards the equator. 

The vertical component {u'zT') is always positive thus 
heat is transported outwards. For A 1 the vertical heat 
fiux increases towards the equator and, in principle, the be- 
havior could be approximated by an expression of the form 
given by equation (|23|l . This is only possible if X|| is negative. 
However, a negative value for X|| is in contradiction with the 
results for the meridional heat fiux which require X|| > 0. 
A positive value for X|| is also obtained from calculations 
of Kitchatinov et al. (1994). Regarding the Coriolis force at 
the pole, n X u — no{~Uy,Ux,0), it is evident that no (lin- 
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Figure 11. Left side: Time and volume averages of the Reynolds stresses Qxz,Qyz and Qzz- Riglit side: Time and volume averages of 
the heat fluxes {u'^T'), {u'yT') and {u'^T'). Both quantities are shown in dependence of the co-latitude 6 for Ta = lO'^. The solid (dotted, 
dashed) curve denotes A = 0.1(1, 10). 



ear) coupling between vertical and horizontal components 
is possible. A transfer of vertical and horizontal momentum 
(through the coupling of Uz with Ux respectively Uy) is only 
possible at higher co-latitude. Therefore one could suppose 
that quantities, like the vertical heat flux or the vertical 
component of the turbulence intensity are maximum at the 
pole and decrease towards the equator. The somehow sur- 
prising opposite behavior shown in the simulations here, is 
confirmed the numerical results of Kapyla et al. (2004) and 



Riidiger et al. (2005). Kapyla et al. (2004) calculated the 
vertical and the meridional heat flux for rotating convection 
in a strongly stratified layer for a various number of rotation 
rates. In all cases they report the maximum of the vertical 
heat flux at the equator and "a tendency" for the minimum 
at a co-latitude 9 « 30°... 45°. These results are confirmed 
by our calculations up to A w 1. Only for strong magnetic 
fields (A « 10) the vertical turbulence intensity and heat 
flux are larger at the pole. 
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4 DISCUSSION AND CONCLUSIONS 

Rotating magnetoconvection has been examined for a wide 
range of parameters and due to the varieties of the results it 
is nearly impossible to describe the outcome uniformly. The 
convectively driven flow changes its characteristic proper- 
ties considerably in dependence of the applied magnetic field 
and/or rotation rate. The behavior of the cell pattern indi- 
cates that in case of a dominant toroidal field convection 
in a fast rotating system should occur in extremely thin 
sheet like structures. In the presented simulations the ver- 
tical extension of the cells is fixed by the vertical boundary 
conditions so that the average size of the convection cells 
in ^-direction is extended over the whole box domain. It is 
hardly believable that such structures are stable over the 
largest possible scales given by the size of the fiuid outer 
core so that the "real" vertical structure of the convection 
remains dubious and might only be obtained from, at present 
unrealistic, global high-resolution simulations. 

For faster rotation the effects of the magnetic field on 
the fluid flow via the Lorentz force become more signifi- 
cant. This has been shown by the behavior of the anisotropy 
functions Ah and Av- Particular attention deserves the ver- 
tical anisotropy. The isotropic behavior on the small scales 
caused by rotation vanishes with the introduction of a mag- 
netic fleld. As a consequence, the vertical component of the 
turbulence intensity dominates the flow which leads to an 
enhanced vertical mixing of the convective layer. Turbulent 
transport in a rotating system is facilitated through the in- 
troduction of a magnetic field. 

The turbulent heat flux has been calculated as an exam- 
ple for a diffusive process that occurs on scales that are not 
resolved in global simulations. Most remarkably, the ver- 
tical heat flux at the pole slightly increases with increas- 
ing magnetic fleld strength (if the rotation is fast enough) 
and remains at a high level even for very strong magnetic 
fields. Thus, cooling of the Earth's interior is enhanced by 
the presence of a dominant toroidal field. Comparison of the 
results from non-rotating magnetoconvection with simula- 
tions that involve rotation and a horizontal (Fig. [Tjl or a ver- 
tically oriented magnetic field (Fig. [9]), respectively, shows 
that an enhancement of the vertical heat flux at high field 
strength only occurs in combination of a horizontally ori- 
ented magnetic field, a sufficient fast rotation rate and close 
to the north-pole. However, the effects of By on the vertical 
heat transport are relatively small compared to the strong 
quenching for A ^ 2{B ^ O.SBeq) in case of an imposed B^. 

Likewise of interest is the behavior of the meridional 
heat ffux which is always negative, indicating that heat is 
transported towards the poles. From comparison with slower 
rotating magnetoconvection we know that the rotational 
quenching of {u'^T') is weaker than for {u'yT') and {u'-,T'), 
so that it might be possible that in a fast rotating system a 
significant amount of heat is transported towards the poles. 
Therfore, the poles should be warmer than the equator. This 
results in a non-radial oriented large-scale temperature gra- 
dient which acts as a source for a meridional large-scale flow. 
Such a meridional flow - if strong enough - interacts with 
the column like large-scale structure of the convective mo- 
tions and might be important for the understanding of the 
dynamo mechanism in the Earth's core (see e.g. Sarson & 
Jones (1999), who identified a fluctuating meridional flow 



as an explanation for the occurrence of dipole reversals) . In- 
deed, global simulations of Olson & Christensen (2002) yield 
a warmer northpole (see their Fig. 5) but it is not clear if this 
effect is a result of turbulent heat transport. A large scale 
meridional flow is also found in the solar convection zone, 
giving rise to important consequences for the solar rotation 
law (Riidiger et al., 2005). 

The distribution of the temperature on the top of the 
fluid core depends on the heat ffow at the core mantle bound- 
ary that is usually prescribed as a boundary condition. Seis- 
mological measurements indicate a non-uniform heat flow at 
the core mantle boundary, whose influence has been exam- 
ined in numerical simulations of Olson & Christensen (2002) . 
The lateral and longitudinal variations of the boundary heat 
ffux affect the local behavior of the ffuid ffow: convection is 
enhanced (inhibited) where the boundary heat ffow is high 
(low). In case of non-axisymmetric boundary conditions at 
the top of the ffuid core, the azimuthal heat ffux becomes 
important. The azimuthal heat ffux is always negative and 
increases towards the equator. This leads to a latitude de- 
pendent westward transport of heat which should interact 
with the constraints given by the inhomogeneous thermal 
conditions at the core mantle boundary. 

Most considerations within the presented work rely on 
the assumption of a dominant toroidal field within the fluid 
core. This does not have to be necessarily the case. Due to 
the insulating mantle the toroidal component of the Earth's 
magnetic field is hidden from any observations. Indeed, most 
of the authors favor a geodynamo-model of aJ?-type (Fearn 
1998) which is characterized by a dominant azimuthal field. 
However, estimations of the toroidal field strength from the 
extrapolation of the dc electric potential near the top of the 
Earth's mantle to the core mantle boundary result in com- 
parable strength of the toroidal and the poloidal field com- 
ponents (Levy & Pearce, 1991). Furthermore, recent simula- 
tions of Christensen & Aubert (2006) that cover a very large 
parameter space, yield solutions without significant differen- 
tial rotation where poloidal and toroidal field strength are 
comparable (a^-Dynamo) . 

Independent from the orientation of the dominant field 
component, the direction of the horizontal components of 
the heat flux (towards the poles and westwards, i.e. against 
the direction of the rotation) is a consequence of the rotation 
and not caused by the magnetic field which only enhances 
or reduces the magnitude of the heat flux. Therefore the 
horizontal components of the heat flux should show a similar 
sign even in case of comparable toroidal and poloidal field 
strength. 

The circumstances within a rotating spherical shell are 
more complicated as they appear in the presented Carte- 
sian model where relevant characteristics have not been con- 
cerned for simplicity. In particular, the combined effects of 
all three components of the turbulent heat ffux cannot be 
inferred from a simple local model so that the consequences 
of the common behavior could only be examined by a global 
simulation. However, the results should at least qualitatively 
resemble the combined effects of the small scale fluctuations 
on large scale flelds and comparable effects are expected to 
be present in spherical geometries. To what extend the re- 
ported results for the heat ffux will endure in simulations in 
spherical geometry, where also possible larger scale ffows will 
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occur, again can only be clarified in highly resolved global 
simulations. 
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